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Abstract 

Nonlinear charge transport parallel to the layers of p-modulation-doped GalnNAs/GaAs quantum wells (QWs) is 
studied both theoretically and experimentally. Experimental results show that at low temperature, T = 13 K, the 
presence of an applied electric field of about 6 kV/cm leads to the heating of the high mobility holes in the 
GalnNAs QWs, and their real-space transfer (RST) into the low-mobility GaAs barriers. This results in a negative 
differential mobility and self-generated oscillatory instabilities in the RST regime. We developed an analytical model 
based upon the coupled nonlinear dynamics of the real-space hole transfer and of the interface potential barrier 
controlled by space-charge in the doped GaAs layer. Our simulation results predict dc bias-dependent self- 
generated current oscillations with frequencies in the high microwave range. 



Introduction 

During the past decade, dilute nitrides, particularly the 
quaternary material system of GalnNAs/GaAs, have 
attracted a great deal of attention, both because of unu- 
sual physical properties and potential applications for a 
variety of optoelectronic devices. The addition of a small 
amount of nitrogen induces a strong perturbation in the 
conduction band of matrix semiconductors, while hav- 
ing a negligible effect on the valence band. As a result, 
the electron mobility is greatly lowered and the hole 
mobility can become higher than the electron mobility, 
in materials with relatively high nitrogen content. High 
hole mobility coupled with the low hole confinement 
energy (110 meV in our calculation for the samples 
investigated in this study) [1] in the GalnNAs/GaAs 
quantum well (QW) structure makes it possible for 
holes in the well to gain enough energy to overcome the 
small band discontinuity under an electric field applied 
parallel to the layer interface, and to transfer into the 
low-mobility /?-doped GaAs layer. This leads to a nega- 
tive differential mobility (NDM) caused by real-space 
hot hole transfer, as we previously observed [1]. There- 
fore, under dc conditions, a self-generated current 



oscillation in the real-space regime, as proposed by 
Scholl and co-authors [2-5], is expected in ^-modula- 
tion-doped GalnNAs/GaAs heterostructures. 

In this work, we study the nonlinear charge transport 
in a modulation-doped GalnNAs/GaAs semiconductor 
heterostructure where the GaAs barrier layer is inten- 
tionally /?-doped. The charge transport processes per- 
pendicular and parallel to the layers far from 
thermodynamic equilibrium are modeled by several 
coupled nonlinear dynamics equations. In this model, 
self-generated current oscillations can be described in 
the following way. Real-space transfer (RST) of holes 
out of the GalnNAs well layer leads to an increase of 
the hole density in the GaAs barrier, which diminishes 
the negative space charge that controls the band bend- 
ing (Figure 1). Consequently, the potential barrier 0 5 
decreases, with some delay due to the finite dielectric 
relaxation time. This leads to an increased thermionic 
emission backward current J b _ w into the GalnNAs well, 
which decreases the hole density in the GaAs barrier. 
As a result, the space charge and ®# are increased in 
the GaAs. This, in turn, decreases the thermionic emis- 
sion backward current from the well into the barrier [6]. 
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Figure 1 Schematic energy-band profile of a GalnNAs/GaAs 
heterostructure. 



Negative differential resistance instabilities in p- 
modulation-doped GalnNAs/GaAs QWs 

The layer structure of the sample used in this study is 
given in Table 1. The sample, which was grown by 
molecular beam epitaxy (MBE) on semi-insulating GaAs 
substrate, consists of three 7 nm thick GalnNAs QWs, 
separated by 20 nm thick Be-doped GaAs barriers. 
These j?-type-doped barriers are separated from the 
QWs by 5 nm undoped spacer layers to reduce the 
remote impurity scattering. The mole fraction of indium 
and nitrogen in the Ga^Jn^NyAs^ QWs is x = 0.3 and 
y - 0.015, respectively. The sample was fabricated in the 
shape of a simple bar for I-V measurement. Fabrication 
details are given somewhere else [1]. 

The nonlinear transport processes depicted in Figure 1 
are modeled by a set of dynamic equations relevant to 
current instabilities in semiconductors. We derive a set 
of nonlinear partial differential equations for the hole 
density in the wells (p w ), and in the barriers (p b ), the 
potential barrier in each of GaAs layers (0#), and the 
dielectric relaxation of the applied parallel field (f//). 

The dynamics of the carrier density in the well and in 
the barrier are given by [5] 



dt 



1 

4 L w 



( J w-b J b-w ) 



(i) 



Table 1 Numerical parameters used in the simulation for 
the GalnNAs/GaAs sample [7] 



Material 



Thickness (A) 



Doping (nv 



GaAs (cap) 


500 


Be: 1 x 10 24 


x3 


GaAs (barrier) 


200 


Be: 1 x 10 24 


x3 


GaAs (spacer) 


50 


UD 


x3 


Ga 1 _ x ln x N y As 1 _ y QW 


70 


UD 


x3 


GaAs (spacer) 


50 


UD 


x3 


GaAs (barrier) 


200 


Be: 1x10 24 


x3 


GaAs (buffer) 


500 


UD 


x3 



9pi L = 
dt 



'{h-w - Jw-b ) 



(2) 



where J w . b and J b _ w are the thermionic currents flow- 
ing from the GalnNAs well layers to the GaAs barrier 
and from the barrier into the well, respectively, q is 
the positive electron charge, and L w and L b are the 
width of the GalnNAs QW and the GaAs barrier, 
respectively. The electric field parallel to the layer 
interface f /7 can be derived from Poisson's equation 
and is given by 



Mil 

dy 



L 1I( + L] 



[{H A -p b )L b +p w L w ] 



(3) 



where s 0 and s s are the absolute and relative permit- 
tivity, respectively. Using Equations (l)-(3), the dielec- 
tric relaxation of the applied parallel field f /7 as a 
function of the current flow (j-direction), the trans- 
verse space coordinator (^-direction), and the time t 
can be written as 



dy 







{ * ) 





d Pu 



dt " dt w 

1 (1 J \ +u gCgtggj 



h+ (4) 
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(Jw-b ~ J b-w ) + Mi 



dx 



L,„ + L h 



Mil 
dy 



Mil 

■ dy 



(5) 



Where <f 0 = il 0 /d is the applied field and il 0 is the 
applied voltage, o L = dl \h(L w +L b )q(A, w R L N A \ is con- 
nected to the load resistance R L , d is the sample 
length, h is the width of the sample, fi w and (A b are the 
hole mobility in the QW and the GaAs barrier, respec- 
tively. By integrating both sides of Equation (5), we 
finally have the dielectric relation of the parallel elec- 
tric field 



i 



Lw + Lb 



Up 



+ WbVb L b]%n (6) 



where /// is the external current density flowing 
through the external circuit at applied bias voltage U 0 . 

Here, we define the current density flowing through 
the sample as a function of applied parallel field, using 



Semi-insulating GaAs substrate 



J n = 



i 



L,„ + Li 



-[qPwVwLw + 4PbVb L b]%H 



(7) 
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The time-dependent potential barrier in the GaAs 
layer is given by 



dt 



e 0 £ s 



4>» + 



(8) 



Equations (1) and (2) represent particle continuity, 
where the thermionic current densities J b . w and J w . b can 
be calculated using Bethe's theory, by assuming that the 
width of the space charge is comparable to the mean 
free path L m of the holes [4,5] 



'w-b 



K B T w 

2nm u 



1/2 



K R T,„ 



(9) 



Table 2 Numerical parameters used in the simulation for 
the GalnNAs/GaAs sample [1] 



Lw 


7 mm 


Lb 


25 mm 


AEv 


0.12 eV 


™w 


0.105 mO 


m h 


0.62 mO 


d 


50 |jm 


h 


28 |jm 


T K 


0.2 ps 




0.1 ps 


TL 


13 K 


UW 


0.3 m2/Vs 


ub 


0.021 m2/Vs 



Jw-b - ~QPb 



K B T b 

2nm h 



1/2 



(10) 



where m w and m h are the hole effective mass in the 

GalnNAs QW and GaAs barrier, respectively, and the 
hole temperatures in the well and barrier are approxi- 
mately given by 

T w~ T L + *eJHwZu T b~ T L + *E b 4Vb%II ( U ) 

Since the number of holes in the well and the barrier are 
related to each other, their total number is conserved [1]: 

Po L w = Pw L w + Pb L b ( 12 ) 

where p 0 is the 3 D hole density in the well at low field. 
Numerical results 

The steady-state can be evaluated by setting Equations 
(1), (2), (6), and (8) to zero, and using the parameters 
listed in Table 2. The resulting static current density 
characteristic as a function of the static electric field is 
shown in Figure 2. The measured 7-1/ curve obtained 
with the same sample in our previous study is placed in 
the figure inset for comparison [1]. Simulation results 
predict that the RST of hot holes leads to an N-shaped 
characteristic with a regime of negative differential resis- 
tance [4,7,8]. The critical field for the onset of NDM is 
the order of 6 kV/cm, which agrees well with our 
experimental results. 

The time-dependent nonlinear Equations (1), (2), (6), 
and (8) have been numerically resolved using Euler's 
methods. The simulation reveals that the instability of 
the dynamic system is strongly dependent on the applied 
dc bias field, <f 0 = U 0 ld. We found that self-generated 
nonlinear oscillation appears in a range of applied dc 



electric fields where the load line lies in the NDM regime, 
as shown in Figure 3a. Figure 3b shows the correspond- 
ing current-density oscillations with frequency of 44 

GHz, for = 10.1 kV/cm and N A = 2.2 x 10 16 cm' 3 . 

It is interesting to find that the oscillation frequency is 
strongly dependent on the dopant concentration in the 
barrier and the barrier thickness, as shown in Figure 4. 
The oscillation frequency increases from 29 to 50 GHz 
as the dopant concentration in the barrier increases 
from 1.9 x 10 16 cm' 3 to 2.4 x 10 16 cm" 3 , accompanied 
by gradually reduced oscillation amplitude. Finally, the 
periodic oscillation damps out when the dopant concen- 
tration is above 2.4 x 10 16 cm" 3 , as shown in Figure 5. 
The oscillation shows similar behavior as the barrier 
thickness increases. The fact that the self-generated 
oscillation frequency can be tuned by the doping con- 
centration and the layer width can be explained by the 
nonlinear combination of the effective thermionic 




0.8 1 1.2 
Electric Field (V/m) 



Figure 2 Static current density-field characteristic as a function 
of the static electric field £ * 7 . The measured l-V characteristic of 
p-modulation-doped sample is shown in the inset. 
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Figure 3 (a) Static current density versus electric field ^curve. The load line (straight line) lies within the NDM area to determine the 
applied dc field, (b) Time-dependent current density curve, with N A = 2.2 x 10 16 cm" 3 . 
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Figure 4 Oscillation frequency as a function of (a) barrier thickness and (b) doping concentration in the GaAs barrier for £ 0 = 24 kV/cm 
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Figure 5 Periodic oscillation damping with N A = 2.4 x 10 16 cm" 3 . 
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emission time, Tq = e 3/2 L w y]37rm* w / AE V and the 

dielectric relaxation time, r r = s 0 s s /^ b A^ A as suggested 
by Dottling and Scholl [9]. The hysteretic switching 
transitions between the stable stationary state and the 
periodic oscillation in a uniform dynamic system depend 
on the ratio of the effective thermionic emission time 
and the dielectric relaxation time, y. In our case, r 0 = 
0.21ps, the change in dopant concentration from 1.8 x 
10 16 cm" 3 to 2.5 x 10 16 cm" 3 leads to y increases from 
0.076 to 0.12 resulting in phase transition in dynamic 
system. 

Conclusion 

In this work, we studied the transport processes parallel 
and perpendicular to the layers of p-type modulation- 
doped GalnNAs/GaAs multi-QW structures far from 
the thermodynamic equilibrium. The simulation results 
of the steady-state predict an NDM induced by RST of 
hot holes in the QWs and the critical electric field of 
the onset of NDM to be the order of 6 kV/cm. This 
value agrees well with our previous experimental results. 
The numerically time-dependent simulations indicate 
that the self-generated oscillation caused by RST with 
the frequency in the range 20-50 GHz appears under 
the right applied electric field. The frequency of self- 
generated oscillation can be flexibly optimized to the 
range of considerable interest for applications as a sim- 
ple way of generating high-frequency microwave power 
based on GalnNAs material system. According to our 
simulation, the predicted self-generated oscillation can 
be observed if the GalnNAs QW structure is optimized 
around 25 nm barrier and less than 2.4 x 10 16 cm" 3 
doping concentration. The current oscillation measure- 
ments will be performed using optimized structures fab- 
ricated into two terminal devices, and shunted with a 50 
H resistor and high-speed circuit (high-speed oscillo- 
scope and pulse generator). The experiment results are 
expected to be published in the near future. 
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